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Abstract 

We present our theoretical results on the mechanism of two transitions in vanadium spinel oxides J4V2O4 {A=Zn, 
Mg, or Cd) in which magnetic V cations constitute a geometrically-frustrated pyrochlore structure. We have derived 
an effective spin-orbital-lattice coupled model in the strong correlation limit of the multiorbital Hubbard model, 
and applied Monte Carlo simulation to the model. The results reveal that the higher-temperature transition is a 
layered antiferro-type orbital ordering accompanied by tetragonal Jahn- Teller distortion, and the lower-temperature 
transition is an antiferromagnetic spin ordering. The orbital order lifts the magnetic frustration partially, and induces 
spatial anisotropy in magnetic exchange interactions. In the intermediate phase, the system can be considered to 
consist of weakly-coupled antiferromagnetic chains lying in the perpendicular planes to the tetragonal distortion. 
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Vanadium spinel oxides AY2O4 with nonmagnetic 
cations A=Zn, Mg, or Cd are one of the most typi- 
cal geometrically- frustrated systems; the magnetic V 
cations constitute the pyrochlore lattice, which is a 
three-dimensional network of corner-sharing tetrahe- 
dra. The geometrical frustration strongly suppresses 
development of magnetic correlations, and the anti- 
ferromagnetic (AF) ordering sets in at a much lower 
temperature ~ 40K than the Curie- Weiss temperature 
^ lOOOK [1]. Besides the AF transition, at a slightly 
higher temperature ~ 50K, the system shows another 
phase transition with the structural change from high- 
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temperature cubic to low-temperature tetragonal sym- 
metry [1] . The origin of the two successive transitions 
is under debate. In particular, theories based on spin 
degree of freedom only or the spin-lattice coupling ap- 
pear to be insufficient to describe the two transitions 
in a comprehensive manner [2] . 

In the present study, we focus on the t2g orbital de- 
gree of freedom in these vanadium spinels to under- 
stand the microscopic mechanism of the two successive 
transitions. Each V^^ cation has two 3d electrons in 
threefold f2g orbitals, and therefore the orbital degree 
of freedom is active. Starting from the standard multi- 
orbital Hubbard model, we consider the perturbation 
in the strong correlation limit to describe the low en- 
ergy physics of these insulating materials, and derive 
the effective spin-orbital-lattice coupled model in the 
form [21 



H = Hso + H]T, 

Hso — —J ^ hij — Js ^ hij, 
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+C(1 — Si • Sj)nia^^j^nj^(^ij^^ (3) 
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where Si is the 5=1 spin operator and riia is the den- 
sity operator for site i and orbital a — 1 (dyz), 2 (dzx), 
3 {da:y)- Here, ftia = 1 — riia and e^ = nn + 71^2 - 2*1^3; 
and we impose a local constraint X]c«=i '^^^ = 2 at 
each site. The summations with {i,j} and {{i,j)) are 
taken over the nearest-neighbor (NN) sites and third- 
neighbor sites, respectively. Here, we take into account 
only the dominant o"-bond hopping integrals in the 
original multiorbital Hubbard model, which results in 
the orbital diagonal interaction in //so; ciiij) is the 
orbital which gives rise to the a bond between sites i 
and j. HjT describes the orbital-lattice coupling part, 
where 7 is the electron-phonon coupling constant of 
the tetragonal Jahn- Teller (JT) mode, Qi denotes the 
amplitude of local lattice distortion at site i, and A 
describes the interaction between NN JT distortions, 
which mimics the cooperative aspect of the JT distor- 
tion. The parameters in Hso are given by the coupling 
constants in the starting multiorbital Hubbard Hamil- 
tonian, and we use the reasonable estimates as J3/J = 
0.02 with J ~ 200K, A = 1.21, B = 0.105, and C = 
0.931. For the JT parameters, we take 7^/ J = 0.04 and 
A/J = 0.15, which are typical values to reproduce the 
tetragonal distortion observed in experiments. 

We have performed Monte Carlo (MC) simulation 
to investigate thermodynamic properties of the model 
(1), which is a classical simulation to avoid the negative 
sign problem due to the geometrical frustration. Since 
quantum nature exists only in the spin operators in 
Hso, we approximate them by classical vectors, and 
apply a standard metropolis MC algorithm. Details of 
MC calculations will be reported elsewhere [3] . 

MC results show that the model (1) exhibits two 
transitions: One is first order transition at To — 0.19J 
corresponding to the orbital ordering with the tetrag- 
onal JT distortion, and the other is at Tn — 0.115J 
with continuous growth of the AF spin ordering. The 
orbital order below To is a layered antiferro type; {d^y, 
dyz) and {dxy, dza;)-occupied planes stack alternatively 
in the c direction (c is the axis of the tetragonal distor- 
tion). This orbital ordering plays a crucial role to re- 
duce the magnetic frustration in the following way: Let 
us consider an effective spin Hamiltonian by replacing 
the orbital part in Hso by its mean-field. The model 
has highly anisotropic spin exchanges; for instance, 
between NN sites, the strong AF interaction Jaf — 
0.931J in the xy direction whereas the weak ferromag- 
netic interaction Jf ~ —0.105 J in the yz and zx direc- 
tions. Moreover, the weak ferromagnetic interactions 
Jp are completely frustrated between the AF xy chains 




Fig. 1. Monte Carlo results of the anisotropy between the 
staggered moments in the xy chains and in the yz ehains. The 
orbital ordering temperature To and the antiferromagnetic 
spin ordering temperature Tn are shown. The system size is 
L X 16 sites. Inset: Anisotropic spin exchange interactions 
under the orbital ordering. 

due to the geometry of the pyrochlore structure. These 
anisotropic exchange interactions are shown in the in- 
set of Fig. 1 in a tetrahedron unit. As a consequence, 
the magnetic correlation becomes highly anisotropic, 
and the system looks like a weakly-coupled AF chains. 
Figure 1 shows the ratio between the staggered mo- 
ments along the xy chains and along the yz{zx) chains. 
Note that the latter moment corresponds to the or- 
dering structure with periodicity four along the yz{zx) 
chains. Above Tn, both moments are zero in the limit 
of 1/ — > 00, and therefore the ratio gives information on 
the anisotropy of the magnetic correlation length. The 
result shows that the anisotropy is suddenly enhanced 
below To and takes a large value for Tn < T < To- 

Below Tn, the AF xy chains are aligned by the third- 
neighbor interaction J3 to form the three-dimensional 
AF order, whose pattern is consistent with the neu- 
tron scattering results [4] . There, the ratio in Fig. 1 ap- 
proaches 1 in the limit of L — > 00 since both moments 
are finite and the same. 

The one-dimensional anisotropy in the intermediate 
phase Tn < T < To can be observed, for example, by 
neutron scattering measurement. Such experiment is 
desired to confirm our scenario. 
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